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ABSTRACT 

We present partial-correlation analyses that examine the strengths of the relationships between Z 2 soo A' h kev, 
a ox , and redshift for optically-selected AGNs. We extend the work of Strateva et al. (2005), that analyzed 
optically-selected AGNs from the Sloan Digital Sky Survey (SDSS), by including 52 moderate-luminosity, 
optically-selected AGNs from the COMBO-17 survey with corresponding deep (w 250 ks to 1 Ms) X-ray ob- 
servations from the Extended Chandra Deep Field-South. The COMBO-17 survey extends ~ 3 magnitudes 
deeper than the SDSS and probes the moderate-luminosity AGNs that numerically dominate the AGN popu- 
lation in the Universe. We also include recently published observations of 19 high-redshift, optically-selected 
AGNs, and 46 luminous, low-redshift AGNs from the Bright Quasar Survey. The full sample used in our 
analysis consists of 333 AGNs, extending out to z ~ 6, with 293 (88%) having X-ray detections. The sample 
spans five decades in UV luminosity and four decades in X-ray luminosity. We confirm that a ox is strongly 
anti-correlated with ^500 A (13.6<r), the highest significance found for this relation to date, and find evidence 
suggesting that the slope of this relation may be dependent on ^500 A- We find that no significant correlation 
exists between a ox and redshift (1.3a), and constrain the maximum evolution of AGN UV-to-X-ray flux ratios 
to be less than 30% (la) out to z — 5. Using our sample's high X-ray detection fraction, we also find a signif- 
icant anti-correlation (3. Oct) between a ox and ^2kev- We make comparisons to earlier studies on this topic and 
discuss implications for X-ray vs. optical luminosity functions. 

Subject headings: Galaxies: Active: Nuclei — Galaxies: Active: Optical/UV/X-ray — Galaxies: Active: 
Evolution — Methods: Statistical 



1. INTRODUCTION 

The relationship between X-ray and UV luminosity 
for AGNs, defined by Tananbaum et al. (1979) as 
a ox = log{l{v x _ n , y )/l(v m )}/log(v x _ T .,Jv m ), where l(v) is the 
monochromatic luminosity (erg s _1 Hz -1 ) at frequency v 
in the rest-frame, has been the subject of many studies over 
the last two decades (e.g., Avni & Tananbaum 1982; Kriss 
& Canizares 1985; Avni & Tananbaum 1986; Anderson & 
Margon 1987; Wilkes et al. 1994; Vignali et al. 2003; Strat- 
eva et al. 2005). Most of these investigations examined 
the relationship between the 2500 A and 2 keV monochro- 
matic luminosities, reducing the a ox relationship to a ox = 
0.3838 log(^2 kev/ ^250oa)- These studies have usually found 
that a ox is anti-correlated with UV luminosity, with at most a 
weak anti-correlation with redshift (but see Yuan et al. 1998; 
Bechtold et al. 2003). Understanding the relationship between 
the intrinsic UV and X-ray emission from AGNs and its evo- 
lution through cosmic time is important for testing energy 
generation models in AGNs, deriving bolometric corrections 
for AGNs which are used in estimating their luminosities and 
accretion rates, identifying X-ray weak AGNs, and estimat- 
ing additional X-ray emission linked with jets in radio-loud 
AGNs. In addition, understanding how a ox evolves with red- 
shift aids in resolving discrepancies between the luminosity 
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functions derived from optical (e.g., Boyle et al. 2000; Wolf 
et al. 2003; Richards et al. 2005) and X-ray (e.g., Miyaji et al. 
2000; Steffen et al. 2003; Ueda et al. 2003; Barger et al. 2005; 
Hasinger et al. 2005; La Franca et al. 2005) AGN samples, 
which map the evolution of the accretion onto supermassive 
black-holes (SMBHs) in the Universe. 

To provide tight constraints upon a ox (l, z), any degenera- 
cies between luminosity and redshift must be broken. This 
is difficult in flux-limited samples, as I and z are typically 
strongly correlated. To obtain good coverage of the I — z 
plane while maintaining a high X-ray detection fraction, it 
is necessary to combine both wide-field and deeper, narrow- 
field surveys. To achieve this goal we expand here upon the 
a ox work of Strateva et al. (2005, hereafter S05), who used the 
wide-field Sloan Digital Sky Survey (SDSS; York et al. 2000) 
quasar sample along with low-redshift (z < 0.2) Seyfert Is 
and high-redshift (z > 4) optically selected AGNs. We add 
to this sample moderate-luminosity, optically-selected AGNs 
from the deep, narrow-field (0.26 deg 2 ) COMBO-17 survey 
(Wolf et al. 2004) centered on the Extended Chandra Deep 
Field-South (E-CDF-S; Lehmer et al. 2005), along with a 
well-studied sample of luminous Bright Quasar Survey (BQS) 
quasars (Schmidt & Green 1983) at z < 0.5 and additional 
optically-selected AGNs at z > 4 (Kelly et al. 2005; Shemmer 
et al. 2005; Vignali et al. 2005). The AGNs optically selected 
in the COMBO-17 survey extend ~ 3 magnitudes fainter than 
those used in previous a ox studies (e.g., B = 19.2; Avni & 
Tananbaum 1986 and B = 19.5; Anderson & Margon 1987) 
and have a substantially higher sky density than the SDSS 
quasar sample (~ 670 deg -2 versus ~ 10 deg -2 ; Wolf et al. 
2004; Schneider et al. 2005). These COMBO-17 AGNs are 
representative of most of the optically-selectable AGNs in the 
Universe. 
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While deep X-ray surveys can provide even higher AGN 
sky densities (up to ~ 7000 deg~ 2 ; e.g., Alexander et al. 
2003), X-ray selected samples are dominated by obscured 
sources (e.g., Bauer et al. 2004). In the UV, this obscura- 
tion leads to increased contamination from the host galaxy 
and uncertain absorption corrections (e.g., Moran, Filippenko, 
& Chornock2002). The X-ray absorption corrections are also 
often uncertain, mainly due to the poor signal-to-noise ratio of 
X-ray spectra with low photon counts. It is therefore difficult 
to measure the intrinsic emission from these obscured AGNs, 
and thus they do not significantly elucidate the intrinsic en- 
ergy generation mechanisms within AGNs. 

The exclusion of obscured sources from our analyses does 
not necessarily mean our conclusions below are valid for un- 
obscured sources only. Our results are still applicable to 
obscured AGNs if the obscuration is a line-of-sight orienta- 
tion effect and does not affect the intrinsic energy generation 
mechanisms within AGNs (i.e., if the "Unified" AGN model 
is valid; see Antonucci 1993, for a review). This will apply 
even given the recent evidence for luminosity-dependent ob- 
scuration found in X-ray selected AGN samples (e.g., Steffen 
et al. 2003; Ueda et al. 2003; Barger et al. 2005; La Franca 
et al. 2005). 

We briefly introduce the optical and X-ray samples used in 
this paper in §2. In §3 we discuss the analysis of the correla- 
tions between ^ 2 soo A' ^2kev, ct ox , and redshift. We discuss our 
findings and present our conclusions in §4. We use J2000 co- 
ordinates and the cosmological parameters Ho = 70 km s _1 
Mpc" 1 , tt M = 0.3, and A = 0.7. 



2. SAMPLE 

To search for correlations among I 



2500 A' 



^2kev, a ox , and red- 
shift within optically-selected AGN samples, it is best to cover 
as much of the I — z plane as possible to minimize degen- 
eracies. It is also important to eliminate AGNs with intrin- 
sic absorption that can affect UV and X-ray photons differ- 
ently, and radio-loud AGNs that can have additional UV and 
X-ray flux associated with the radio jet (e.g., Worrall et al. 
1987; Wilkes & Elvis 1987; Worrall & Birkinshaw 2005). To 
this end, we combine the data from S05 with the optically- 
selected, moderate-luminosity AGNs from the COMBO- 17 
survey, a subsample of sources from the BQS quasar catalog 
which includes higher-luminosity, optically-selected AGNs, 
and additional optically-selected, z > 4 AGNs. We describe 
our AGN sample below starting with a brief description of the 
S05 catalog. 

2.1. Strateva et al. (2005) Sample 

The main S05 sample is composed of broad-line AGNs 
optically selected via multi-band photometry in the SDSS 
Data Release 2 (SDSS DR2; Abazajian et al. 2004), exclud- 
ing broad absorption-line quasars (BALQSOs) and radio-loud 
sources. The S05 sample is an unbiased subsample of the 
full SDSS DR2 AGN catalog, selecting sources with medium- 
deep (> 11 ks) ROSAT PSPC observations. Of the 155 SDSS 
DR2 AGNs in the S05 sample, 126 (81%) have X-ray detec- 
tions. X-ray upper limits were calculated for the remaining 29 
sources. The 155 AGNs span redshifts 0.2 < z < 3.5. 

To extend their I — z plane coverage, S05 include the low- 
redshift Seyfert 1 sample of Walter & Fink (1993), excluding 
the radio-loud sources (Lsghz > 10 25 W Hz -1 ). Twelve of 
these AGNs are also included in the well-defined BQS (see 
§ 2.2.2), and thus we have removed them from the S05 sam- 
ple. For their high-redshift sample, S05 include optically- 
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FIG. 1 . — Apparent i-band PSF magnitude vs. relative g — i color, A(g — 
i), for our selected COMBO-17 sample (red diamonds), the SDSS sample 
(blue circles), and the BQS sample ( vioiet downward-pointing triangles). For 
comparison the full SDSS DR2 sample is shown as a blue contour enclosing 
90% of the data and small blue dots representing the outliers. Any sources 
with significant galaxy contributions that could affect their observed colors 
(identified spectroscopically for SDSS AGNs and labeled 'QSO (GAL?)' in 
the COMBO-17 survey) are enclosed by large squares. COMBO-17 sources 
with z > 2.1 are denoted by open diamonds. The BQS sample sources 
observed in the SDSS are denoted with large symbols, and sources with SDSS 
magnitudes derived from other bands are denoted with small symbols. In the 
computation of A(g — i), only AGNs with point-source morphology were 
used to determine the median g — i color as a function of redshift to prevent 
artificial reddening due to host-galaxy contamination. The characteristic i- 
band magnitude ranges for the BQS, SDSS, and COMBO-17 samples are 
denoted along the left-hand side of the figure. 

selected, z > 4 quasars from the SDSS, the Palomar Digital 
Sky Survey (Djorgovski et al. 1998), and the Automatic Plate 
Measuring facility survey (Irwin et al. 1991). These high- 
redshift AGNs were specifically targeted by Chandra and 
XMM-Newton, and S05 again eliminated radio-loud sources 
and BALQSOs from their sample. It is important to note 
that all of these optically-selected, high-redshift AGNs would 
have met the SDSS AGN color-selection criteria, so the intro- 
duction of the sources from other surveys does not bias the 
S05 sample. This high-redshift sample includes 36 AGNs, in- 
creasing the S05 total sample size to 216 optically-selected 
AGNs (after removing the twelve Seyfert Is also observed in 
the BQS), 183 (85%) of which have detected X-ray counter- 
parts and 33 (15%) of which have X-ray upper limits. 

2.2. New Objects in This Study 

2.2.1. COMBO-17 

To extend the AGN sample to fainter luminosities we in- 
clude the optically-selected AGNs from the portion of the 
COMBO-17 survey that covers the E-CDF-S (Wolf et al. 
2004). To determine if the COMBO-17 AGNs are indeed con- 
sistent with being a fainter extension of the SDSS AGNs, we 
examine the color-magnitude distributions of the two popu- 
lations. Figure 1 shows the apparent i-band PSF magnitudes 
versus the relative g — i colors for the COMBO-17 AGNs 
(diamonds), the SDSS AGNs (circles), and the BQS sample 
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covered by SDSS (large triangles) and with SDSS colors de- 
rived from APM magnitudes (small triangles). The SDSS g 
and i magnitudes were calculated for the COMBO- 17 AGNs 
using the observed BVRI magnitudes and the color transfor- 
mations calculated by Jester et al. (2005). These AGN color 
transformations are for sources with redshifts z < 2.1. In Fig- 
ure 1 we denote COMBO- 1 7 sources beyond this redshift with 
open diamonds. The relative g — i colors, denoted as A(g — i), 
were constructed by subtracting the median g — i color of the 
DR2 AGNs as a function of redshift from the observed g — i 
color. Relative (i.e., redshift-corrected) colors can be used to 
detect significant differences in the optical/UV continuua of 
AGNs (Richards et al. 2003), so they can be used to compare 
the SDSS and COMBO-17 AGN populations. From Figure 1 
it appears that the COMBO-17 AGN sample is plausibly the 
faint extension of the SDSS AGN population. 

The COMBO-17 survey used observations in 5 broad-band 
and 12 medium-band filters spanning the optical regime to de- 
rive photometric redshifts and "fuzzy" spectra for over 63 000 
sources in the E-CDF-S field. Wolf et al. (2004) find 175 
sources with photometric colors that are best matched by the 
SDSS broad-line QSO template spectrum of Vanden Berk 
et al. (2001). Their method preferentially selects broad-line 
AGNs, but narrow-line sources with strong AGN emission 
lines could also be selected. 

The 17-band photometry and use of multiple galaxy/ AGN 
templates allows the COMBO-17 survey to measure among 
the most accurate photometric redshifts to date. To examine 
the quality of the photometric redshifts for the COMBO-17 
AGNs, we compared the photometric redshifts with a compi- 
lation of publicly available spectroscopic redshifts of sources 
within the GOODS-South field created by Alessandro Ret- 
tura. 6 For our study, we find spectroscopic redshifts for 12 
COMBO-17 AGNs and use photometric redshifts for the re- 
maining sources. Overall, we find excellent agreement be- 
tween the photometric and spectroscopic redshifts for AGNs 
at the magnitudes we consider below. The spectroscopic red- 
shifts for 1 1 sources differ by less than 5% from their mea- 
sured photometric redshifts. Only one of the twelve COMBO- 
17 sources with measured spectroscopic redshifts has a pho- 
tometric redshift that significantly disagrees with the spec- 
troscopic value (COMBO-17 ID: 30792; z phot = 1.929; 
z spcc = 0.743). 

To obtain monochromatic 2 keV X-ray luminosities (or up- 
per limits) for these COMBO-17 AGNs, we match the sources 
to the E-CDF-S and the Chandra Deep Field-South (CDF-S; 
Alexander et al. 2003) X-ray catalogs. The E-CDF-S is com- 
posed of four 250 ks Chandra observations that cover the en- 
tire COMBO-17 field (Lehmer et al. 2005). The 1 Ms CDF-S 
probes fainter fluxes than the E-CDF-S, but it does not cover 
the entire COMBO-17 field. We use a 2" matching radius 
to identify the X-ray counterparts of the COMBO-17 AGNs, 
and we estimate the false-match probability to be < 0.5%. 
All of the COMBO-17 sources with X-ray counterparts were 
detected in the E-CDF-S, and none had more than one X-ray 
counterpart. We use the CDF-S to derive X-ray upper limits 
only when COMBO-17 AGNs are not detected in the overlap- 
ping E-CDF-S fields. 

To minimize potential contamination of our correlation 
analysis below we impose two unbiased selection criteria on 
the full COMBO-17 AGN sample. First, we only consider 
sources brighter than R = 23. At these magnitudes the prob- 
ability of a COMBO-17 source being misidentified as an AGN 
is low (~ 1%; Wolf et al. 2004). In addition, at R > 23 the 
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FIG. 2.— Histogram of the ctox values for the 52 COMBO-17 AGNs 
included in our sample. The mean of the distribution is shown with a 
dotted line. COMBO-17 AGNs with T < 1.6 are marked (hatched his- 
togram). Sources with a x upper-limits (i.e., without X-ray detections) are 
also marked (shaded histogram). The left-pointing arrow indicates the direc- 
tion of the limits. For comparison, we show the mean q x values for the three 
AGN samples in S05, which are described in Table 3. 

fraction of COMBO-17 AGNs with X-ray detections drops 
quickly. This magnitude threshold ensures that the fraction 
of AGNs with X-ray detections remains high, reducing the 
impact of X-ray limits in our analysis. Second, we only in- 
clude sources that have HST ACS data from either the Galaxy 
Evolution from Morphology and SEDs survey (GEMS; Rix 
et al. 2004) or the Great Observatories Origins Deep Survey 
(GOODS; Giavalisco et al. 2004). These high-resolution data, 
which cover ~ 84% of the COMBO-17 field, can be used 
to minimize the galaxy contribution to the UV luminosity of 
the COMBO-17 sources identified as being extended objects 
(flagged as 'QSO (GAL?)' in Wolf et al. 2004). These two se- 
lection criteria reduce the COMBO-17 AGN sample from 175 
to 60 sources, with 49 (82%) having X-ray detections and 13 
(22%) identified as extended objects. 

Using the high-resolution GEMS and GOODS ACS 
U-band images we identified four COMBO-17 AGNs 
(COMBO-17 IDs: 5498, 18324, 31898, and 50415) that 
were incorrectly classified as extended sources but are actu- 
ally close sources blended together in the lower-resolution 
COMBO-17 images. We correct the COMBO-17 classifica- 
tion for these blended sources. For the remaining nine ex- 
tended COMBO-17 AGNs, we attempted to minimize the 
contribution of the host galaxies to the measured optical/UV 
luminosities of the AGNs. Using the high-resolution ACS im- 
ages we compared the magnitudes of the COMBO-17 AGNs 
using 3" and 0.5" apertures. A mean magnitude difference 
of 0.352 mag was calculated for the unresolved COMBO-17 
AGNs. For the extended sources, we calculated the magnitude 
differences using the same aperture sizes. Any magnitude off- 
sets exceeding the mean difference calculated for the point- 
like AGNs were considered to be caused by excess light from 
the host galaxy. We subtracted this host-galaxy contribution 
from the observed COMBO-17 i?-band magnitude and found 
that six of the nine extended sources subsequently fell below 
our R = 23 magnitude limit. Since these extended AGNs 
would not have met our magnitude criterion without the ad- 

6 Available at http://www.eso.org/science/goods/ 
spectroscopy/CDFS_Mastercat/ 
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TABLE 1 
COMBO-17/E-CDF-S AGN DATA. 



COMBO-17 E-CDF-S a.72000 6.72000 2 log ^2500A- ) log ^2500A- ) lo s(/2kev) lo s( ; 2kev) "ox T R 

ID ID mtcm -3 » -1 Mt -1 l [a»» -1 rn»an -2 » _l «t- l l ro»« -l Rt -l l 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 
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33644 


525 


03 


32 59.85 


-27 47 48.2 


2. 


.538 


-28.15 


30 


.02 




-31 


.74 


26 


76 




-1 


.25 


1 


.71 


< 2.35 


34357 


304 


03 


32 08.67 


-27 47 34.2 





.543° 


-27.33 


29 


.55 




-30 


.55 


25. 


.89 




-1 


.40 


1 


.73 


13.89 


36120 


180 


03 


31 49.41 


-27 46 34.0 


2. 


.306 


-28.35 


29 


.74 




-32 


.09 


25. 


.85 




-1 


.49 


1 


.42 


< 5.10 


37487 


422 


03 


32 39.08 


-27 46 01.8 


1. 


.216° 


-27.86 


29 


73 




-32 


.68 


25 


01 




-1 


.81 


1 


.07 


< 1.25 


38551 


390 


03 


32 29.98 


-27 45 29.8 


1. 


221° 


-28.03 


29 


.56 




-31 


.53 


26 


00 




-1 


.37 


1 


.66 


< 1.88 


38905 


549 


03 


33 03.62 


-27 45 18.7 


1. 


264 


-28.71 


28. 


90 




-32 


.05 


25 


66 




-1 


.25 


1 


.55 


<11.94 


39432 


393 


03 


32 30.22 


-27 45 04.5 


0. 


.735° 


-28.52 


28. 


.63 




-31 


.85 


25 


33 




-1 


.27 


1 


.88 


< 4.83 


42601" 


515 


03 


32 59.07 


-27 43 39.5 


0. 


.510 


-27.98 


28. 


85 




-31. 


.52 


24 


90 




-1 


.51 


1 


.67 


< 3.60 


43151 


249 


03 


32 00.36 


-27 43 19.5 


1. 


.037° 


-28.42 


29 


.03 




-31 


.70 


25 


.48 




-1 


.36 


1 


.54 


< 4.70 


44126 


583 


03 


33 08.78 


-27 42 54.4 





.729 


-28.91 


28 


.23 




-31 


.78 


25. 


.06 




-1 


.22 


1 


.74 


< 11.20 


47615 


191 


03 


31 50.96 


-27 41 15.7 





.649 


-28.52 


28 


.52 




-31 


.96 


25 


.11 




-1 


.31 


1 


.79 


< 5.92 


48284 


378 


03 


32 27.01 


-27 41 04.9 





.734° 


-27.16 


29 


.99 




-30 


.61 


26 


25 




-1 


.43 


1 


.79 


5.59 


48870 


712 


03 


33 28.94 


-27 40 43.7 


2. 


.146 


-28.46 


29 


57 




-31 


.68 


26 


.66 




-1 


.12 


1 


.61 


< 5.64 


49298 


374 


03 


32 26.49 


-27 40 35.5 


1. 


.031° 


-27.52 


29 


93 




-31 


37 


26 


.03 




-1 


.49 


2 


.08 


< 0.61 


50415" 


412 


03 


32 37.45 


-27 40 00.1 


0. 


.666' 


-28.89 


28. 


.17 




-32 


.21 


24 


73 




-1 


.32 


1 


.83 


<11.25 


50997 


416 


03 


32 38.12 


-27 39 44.8 


0. 


.837° 


-27.71 


29 


.55 




-31 


.38 


25. 


.90 




-1 


.40 


2 


.05 


< 0.92 


51835 


118 


03 


31 40.05 


-27 39 17.6 


2. 


.179 


-28.60 


29 


.45 




-32 


.51 


25 


.63 




-1 


.47 


1 


.50 


< 7.73 


52963 




03 


32 52.60 


-27 38 46.2 





.548 


-28.67 


28. 


.22 


< 


-33 


.30 


< 23 


.59 


< 


-1 


.78 






< 7.38 


54839 


116 


03 


31 39.25 


-27 37 52.0 


1. 


.428 


-28.34 


29 


.38 




-31 


.88 


26 


.04 




-1 


.28 


1 


.88 


< 4.61 


54969" 




03 


32 49.28 


-27 37 56.7 


3 


.361 


-27.39 


30 


97 


< 


-33 


09 


< 25 


.61 


< 


-2 


.06 






< 1.23 


56074 


320 


03 


32 11.65 


-27 37 25.9 


1. 


.574 


-26.96 


30 


.83 




-31 


.13 


26 


.84 




-1 


.53 


1 


.44 


19.50 


57653 


35 


03 


31 23.53 


-27 36 31.6 


1. 


.653 


-28.25 


29 


.58 




-32 


.70 


25. 


.25 




-1 


.66 


1 


.40 


< 3.74 


58478 


516 


03 


32 59.19 


-27 36 11.7 


1. 


.348 


-28.01 


29 


.66 




-32 


.05 


25. 


.77 




-1 


.49 


2 


.12 


< 2.55 


60939 


245 


03 


31 58.91 


-27 35 16.1 


2. 


794 


-28.29 


29 


.94 




-32 


.37 


26 


.21 




-1 


.43 





.96 


< 7.46 



a Source is blended with nearby sources in COMBO-17 images. 
b Source is identified as 'QSO (GAL?)' in the COMBO-17 survey. 
c Spectroscopic redshift 

ditional flux from the host galaxy, we exclude them from our 
sample, reducing it to 54 sources. None of these six AGNs 
has an X-ray detection. 

We attempted to remove both radio-loud and obscured 
AGNs from our sample of COMBO-17 AGNs. Using 
1.4 GHz VLA radio maps covering the COMBO-17 field 
(Kellermann et al. 2004) we found that only five of the 54 
COMBO-17 AGNs have detected radio counterparts. We 
used the 1.4 GHz flux limit of 40/iJy (5a), quoted by Keller- 
mann et al. (2004), as an upper limit for the undetected ra- 
dio sources. We calculated the radio-loudness parameter, 
R, using the COMBO-17 i?-band monochromatic flux (A = 
4350 A) and converting the observed 1 .4 GHz flux density to 
rest-frame 5 GHz, assuming S(i') oc ^~ - 8 . We eliminate the 



two COMBO-17 AGNs with R > 30, and give the values 
(or limits) of R for the remaining sources in Table 1 . Our fi- 
nal COMBO-17 AGN sample consists of 52 sources, with 47 
(90%) having X-ray counterparts. 

We interpolate the optical data to find the monochromatic 
luminosity at 2500 A for each COMBO-17 source using the 
redshifts and the 12 medium-band photon fluxes given by 
Wolf et al. (2004). These bands range between 4180 - 9140 A 
which, for rest-frame 2500 A, corresponds to a redshift range 
of 0.67 < z < 2.66. For the 15 sources outside of this red- 
shift range we linearly extrapolate to find the 2500 A lumi- 
nosity. We calculated the monochromatic 2 keV luminosi- 
ties from the 0.5 — 2 keV luminosities, assuming a power-law 
spectral energy distribution (SED) with T = 2.0. Table 1 
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TABLE 2 
BQS/ROSATAGN DATA. 



Name 


Ct T 9(1(1(1 


79(1(1(1 


z 


v J 2500 A/ 

[erg cm -2 s — 1 Hz -1 ] 


&V"2500 A' 

r er9 ^— 1 Hz - 1 1 
L cl & s ni J 


logf Vpv) 
& V J L Kc V / 

[ergcm~ 2 s~ 1 Hz~ 1 ] 


\os.d'} Vf»v) 
"-'feV "Z Ke V / 

L cl & & ni J 














(1) 


(2) 


(31 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


PG 0026+129 


00 29 13.8 


+13 16 05 


0.142 


-25.47 


30.20 


-29.21 


26.46 


-1.44 


PG 0050+124 


00 53 34.9 


+12 41 36 


0.061 


-25.24 


29.69 


-28.96 


25.96 


-1.43 


PG 0052+251 


00 54 52.2 


+25 25 39 


0.155 


-25.51 


30.24 


-29.12 


26.64 


-1.39 


PG 0157+001 


01 59 50.25 


+00 23 40.8 


0.164 


-25.73 


30.07 


-29.83 


25.97 


-1.57 


PG 0804+761 


08 10 58.5 


+76 02 43 


0.100 


-25.01 


30.35 


-28.91 


26.45 


-1.50 


PG 0838+770 


08 44 45.3 


+76 53 10 


0.131 


-25.79 


29.81 


-29.76 


25.84 


-1.52 


PG 0844+349 


08 47 42.47 


+34 45 04.4 


0.064 


-25.40 


29.57 


-30.28 


24.69 


-1.87 


PG 0923+201 


09 25 54.7 


+19 54 04 


0.190 


-25.64 


30.29 


-29.59 


26.34 


-1.52 


PG 0947+396 


09 50 48.39 


+39 26 50.5 


0.206 


-26.05 


29.96 


-29.60 


26.41 


-1.36 


PG 0953+414 


09 56 52.39 


+41 15 22.2 


0.239 


-25.50 


30.64 


-29.50 


26.64 


-1.54 


PG 1012+008 


10 14 54.90 


+00 33 37.4 


0.185 


-25.83 


30.08 


-30.02 


25.89 


-1.61 


PG 1048+342 


10 51 43.8 


+33 59 26 


0.167 


-25.94 


29.88 


-29.85 


25.97 


-1.50 


PG 1049-005 


10 51 51.50 


-00 51 17.7 


0.357 


-25.76 


30.74 


-30.18 


26.32 


-1.70 


PG 1114+445 


11 17 06.40 


+44 13 33.3 


0.144 


-25.88 


29.81 


-30.14 


25.54 


-1.64 


PG 1115+407 


11 18 30.29 


+40 25 54.0 


0.154 


-26.06 


29.69 


-29.86 


25.89 


-1.46 


PG 1116+215 


11 19 08.8 


+21 19 18 


0.177 


-25.26 


30.62 


-29.35 


26.52 


-1.57 


PG 1121+422 


11 24 39.18 


+42 01 45.0 


0.234 


-26.15 


29.97 


-29.88 


26.24 


-1.43 


PG 1151+117 


11 53 49.27 


+11 28 30.4 


0.176 


-25.90 


29.96 


-29.56 


26.30 


-1.40 


PG 1202+281 


12 04 42.1 


+27 54 12 


0.165 


-26.07 


29.73 


-29.52 


26.29 


-1.32 


PG 1211+143 


12 14 17.7 


+14 03 13 


0.085 


-25.04 


30.18 


-28.85 


26.37 


-1.46 


PG 1216+069 


12 19 20.93 


+06 38 38.5 


0.334 


-25.72 


30.72 


-29.58 


26.86 


-1.48 


PG 1229+204 


12 32 03.6 


+20 09 30 


0.064 


-25.48 


29.49 


-29.15 


25.82 


-1.41 


PG 1259+593 


13 01 12.93 


+59 02 06.8 


0.472 


-25.69 


31.06 


< -30.20 


< 26.55 


< -1.73 


PG 1307+085 


13 09 47.00 


+08 19 48.2 


0.155 


-25.54 


30.21 


-29.47 


26.29 


-1.51 


PG 1322+659 


13 23 49.52 


+65 41 48.2 


0.168 


-25.83 


29.99 


-29.56 


26.26 


-1.43 


PG 1352+183 


13 54 35.6 


+18 05 18 


0.158 


-25.78 


29.99 


-29.59 


26.18 


-1.46 


PG 1354+213 


13 56 32.8 


+21 03 51 


0.300 


-26.29 


30.06 


-29.98 


26.37 


-1.42 


PG 1402+261 


14 05 16.2 


+25 55 34 


0.164 


-25.62 


30.19 


-29.45 


26.35 


-1.47 


PG 1404+226 


14 06 21.9 


+22 23 47 


0.098 


-26.10 


29.25 


-29.51 


25.83 


-1.31 


PG 1415+451 


14 17 00.82 


+44 56 06.4 


0.114 


-25.95 


29.52 


-29.97 


25.51 


-1.54 


PG 1416-129 


14 19 03.8 


-13 10 45 


0.129 


-25.12 


30.47 


-29.17 


26.42 


-1.55 


PG 1426+015 


14 29 06.57 


+01 17 06.1 


0.086 


-25.13 


30.10 


-28.91 


26.32 


-1.45 


PG 1427+480 


14 29 43.07 


+47 47 26.2 


0.221 


-26.09 


29.98 


-29.96 


26.11 


-1.49 


PG 1435-067 


14 38 16.2 


-06 58 20 


0.129 


-25.39 


30.20 


-29.47 


26.12 


-1.56 


PG 1440+356 


14 42 07.47 


+35 26 23.0 


0.077 


-25.52 


29.61 


-29.03 


26.10 


-1.35 


PG 1444+407 


14 46 45.94 


+40 35 05.8 


0.267 


-25.85 


30.39 


-29.98 


26.26 


-1.59 


PG 1519+226 


15 21 14.3 


+22 27 44 


0.137 


-25.86 


29.78 


-29.95 


25.69 


-1.57 


PG 1543+489 


15 45 30.24 


+48 46 09.1 


0.400 


-26.00 


30.61 


-30.47 


26.14 


-1.72 


PG 1552+085 


15 54 44.58 


+08 22 21.5 


0.119 


-25.81 


29.71 


-30.13 


25.38 


-1.66 


PG 1612+261 


16 14 13.20 


+26 04 16.2 


0.131 


-25.77 


29.83 


-29.15 


26.45 


-1.30 


PG 1613+658 


16 13 57.2 


+65 43 10 


0.129 


-25.78 


29.81 


-29.09 


26.50 


-1.27 


PG 1617+175 


16 20 11.2 


+17 24 28 


0.114 


-25.32 


30.16 


-29.67 


25.81 


-1.67 


PG 1626+554 


16 27 56.2 


+55 22 32 


0.133 


-25.87 


29.74 


-29.34 


26.27 


-1.33 


PG 2130+099 


21 32 27.82 


+10 08 19.2 


0.061 


-25.20 


29.72 


-28.96 


25.97 


-1.44 


PG 2214+139 


22 17 12.26 


+14 14 20.9 


0.067 


-25.44 


29.57 


-30.66 


24.35 


-2.00 


PG 2233+134 


22 36 07.68 


+13 43 55.3 


0.325 


-25.90 


30.52 


-30.03 


26.39 


-1.59 



lists the properties of the COMBO- 17 AGN sample, includ- 
ing the COMBO-17 ID from Wolf et al. (2004), the E-CDF-S 
ID from Lehmer et al. (2005), the sources' positions and red- 
shifts, monochromatic 2500 A and 2 keV fluxes and luminosi- 
ties, a ox , the effective X-ray power-law photon index (r), R, 
and a column identifying blended or extended COMBO-17 
AGNs. 

Since the majority of the COMBO-17 AGNs do not have 
optical spectra, we cannot easily identify obscured AGNs (us- 
ing emission-line and absorption-line diagnostics) that may 
contaminate our sample. To identify potential obscured 
COMBO-17 AGNs, we will test the a ox correlations below 
after removing the 22 sources with an effective X-ray power- 
law photon index T < 1.6, which is indicative of X-ray 
obscuration (Lehmer et al. 2005). This procedure does not 
help to identify potential obscured AGNs in the small num- 
ber of COMBO-17 sources without X-ray counterparts; these 
sources will need future spectroscopic observations to iden- 
tify possible absorption. 

The a ox distribution of our COMBO-17 AGN sample is 



shown in Figure 2. COMBO-17 AGNs that have measured 
effective photon indices of T < 1.6, and sources that have a ox 
upper-limits, are shown as hatched and shaded histograms, re- 
spectively. The COMBO-17 AGNs typically have flatter a ox 
values than those calculated for more luminous AGN sam- 
ples. This is consistent with the observed anti-correlation be- 
tween a ox and optical/UV luminosity seen in the studies men- 
tioned in Section 1, where the value of a ox falls with increas- 
ing Z 2 5oo A- F° r comparison, we mark the mean a ox values for 
the three AGN samples used in S05 in Figure 2. 

2.2.2. Bright Quasar Survey 

To increase the luminosity range covered at low redshifts 
we included a set of 46 BQS quasars with z < 0.5 ana- 
lyzed by Brandt, Laor, & Wills (2000). While the BQS sam- 
ple does extend beyond z = 0.5, Jester et al. (2005) found a 
bias against detecting AGNs in the BQS catalog in the redshift 
range 0.5 < z < 1.0. Since the BQS sample is meant to be a 
higher-luminosity complement to the S05 and Seyfert 1 sam- 
ples at low redshifts, we removed 21 AGNs with M B > -23. 
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This luminosity threshold also removes fainter AGNs that are 
more likely to be affected by contamination from their host 
galaxies. We excluded 13 radio-loud AGNs and five known 
BALQSOs (see footnote 4 in Brandt et al. 2000). In ad- 
dition, following the arguments of Brandt et al. (2000), we 
also exclude two BQS quasars that have significant UV ab- 
sorption, but do not meet the formal criteria for BALQSOs: 
PG 1351 + 640 (e.g., Zheng et al. 2001) and PG 1411 + 442 
(e.g., Malkan et al. 1987). 

We converted the / 3000 a. values used in Brandt et al. (2000) 
to £2500 A' assuming f (u) oc v~ a with a = —0.5. We used 
PIMMS 7 to calculate the monochromatic 2 keV luminos- 
ity, assuming a power-law spectrum with T = 2, from the 
ROSAT 0.5 - 2.0 keV pointed PSPC count-rate when avail- 
able, and from the 0.5 - 2.0 keV ROSAT All Sky Survey 
(RASS) count-rate when no pointed ROSAT observation ex- 
isted. An X-ray upper limit for one BQS quasar with no X-ray 
counterpart was calculated from the RASS. Our BQS sam- 
ple contains 46 sources, with 45 (98%) having X-ray detec- 
tions. We use the astrometry presented in Jester et al. (2005) 
for the BQS sources included therein and the astrometry from 
Veron-Cetty & Veron (2003) for the remaining sources. Ta- 
ble 2 presents the BQS sources' names, positions, and red- 
shifts along with the monochromatic 2500 A and 2 keV fluxes 
and luminosities, and a ox - 



2.2.3. Additional High-z AGNs 

To increase the size of our high-redshift sample we included 
recently published, optically-selected AGNs with z > 4. We 
included the ten luminous quasars listed in Table 3 of Vig- 
nali et al. (2005) (excluding the BALQSO PSS 1506 + 5220), 
and three radio-quiet AGNs from Shemmer et al. (2005) that 
were not included in S05 (Q 0000 - 263, BR 0351 - 1034, 
and BR 2237 — 0607). To extend the coverage of our sam- 
ple to somewhat lower luminosities at z > 4 we included 
a sample of six optically-selected, faint, radio-quiet, non- 
BALQSO sources targeted in > 10 ks Chandra observations 
(Kelly et al. 2005). We analyzed these X-ray and optical 
data in the same manner as outlined in Vignali et al. (2005). 
We briefly present these data in the Appendix. The addition 
of these three optically-selected, z > 4 AGN samples adds 
19 optically-selected AGNs to our total sample, 18 (94%) of 
which have X-ray detections. 

For comparison purposes, we also plot the eight 8 published 
X-ray-selected, radio-quiet, z > 4, AGNs (crosses) in Fig- 
ures 3—9. We present the observed properties of these sources 
in Table A2 in the Appendix. While these X-ray-selected 
AGNs substantially extend the coverage of the UV luminosity 
- redshift plane at high redshifts (see Figure 3a), we do not in- 
clude these AGNs in our formal analyses as they have been se- 
lected with different criteria than the other optically-selected 
AGN samples. However, in §4.3 we discuss how these AGNs 
fit into the relations derived in §3 using the optically-selected 
AGN samples. 

7 Portable Interactive Multi-Mission Simulator available at 

http://heasarc.gsfc.nasa. gov/ Tool s/w3pimms . html 

8 We excluded the two known X-ray-selected, radio-loud, z > 4 AGNs, 
RX J1028.6-0844 (z = 4.28; Zickgraf et al. 1997) and RX J1759.4+6638 
(z = 4.32; Henry et al. 1994). An up-to-date compilation of z > 4 AGNs 
with X-ray detections is available at http://www.astro.psu.edu/ 
user s/niel /paper s/highz-xr ay-detected . dat 



TABLE 3 

Summary of Samples Utilized 



Sample 


Total 


Number 0.5 


- 2.0 keV Limit 


Area 




AGNs 


X-ray Detected [erg cm 2 s : ] 


[sq. deg.] 




S05 


SDSS 


155 


126 (81%) 


~ io- 14 


15 


Seyfert 1 


25 a 


25 (100%) 


~ io- 12 


b 


High-2 


36 


32 (89%) 


~ io- 15 


. . . b 


This work 


COMBO-17 


52 


47 (90%) 


~ IO" 16 


0.26 


BQS 


46 


45 (98%) 


~ io- 13 


10,714 


High-,z 


19 


18 (95%) 


~ io- 15 


. . . b 


Full Sample 


Total 


333 


293 (88%) 




b 



"After removing the twelve sources that overlap with the BQS sample. 
b Not well defined. 



2.3. Full Sample 

The full sample used in our analysis consists of 333 AGNs, 
293 (88%) of which have X-ray detections; this is a much 
higher X-ray detection fraction than for most previous analy- 
ses, where the detection fraction is typically 10 — 50%. The 
number of AGNs contributed by each sample, and their re- 
spective X-ray detection fractions, are summarized in Table 3. 
Figure 3 shows the monochromatic 2500 A (left) and 2 keV 
(right) luminosity versus redshift for the combined AGN sam- 
ple. The addition of the COMBO-17 (diamonds) and BQS 
(downward-pointing triangles) AGNs to the S05 data substan- 
tially improves the luminosity range covered out to z ~ 3, 
corresponding to ~ 85% of cosmic history. 

3. PARTIAL CORRELATION AND LINEAR REGRESSION ANALYSES 

To measure properly the correlations among Z 2 5oo A' '2 keV, 
a ox , and redshift, we must attempt to eliminate (or minimize) 
potential sources of bias from our analysis. Following the ar- 
guments of Kembhavi, Feigelson, & Singh (1986), we use the 
AGN luminosities in our calculations and not the measured 
fluxes which, if used, can yield different results depending 
upon the cosmology and fc-corrections assumed. However, 
using the luminosities can introduce a different bias since lu- 
minosity and redshift are typically correlated in flux-limited 
samples. Our sample was constructed to cover a large area of 
the / — z plane, which acts to break this degeneracy. 

While the majority of the AGNs in our sample have X-ray 
detections (88%), we cannot overlook the effects of sources 
with X-ray upper limits. In addition, our sample data un- 
doubtedly contain intrinsic scatter that is not quantified in 
our measurements. We have minimized the effects of some 
sources of intrinsic scatter such as host-galaxy contamination, 
obscuration, and excess UV and X-ray flux associated with ra- 
dio jets (as discussed in §2.2.1). Another source of intrinsic 
scatter in our data is caused by AGN variability, which affects 
our calculations of a ox since the optical and X-ray observa- 
tions were not taken concurrently. Unfortunately, the effects 
of variability cannot be corrected and remain a significant 
source of intrinsic scatter in our data. To measure the afore- 
mentioned correlations among these AGNs, proper statistical 
tools must be used that take these issues into consideration. 
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FIG. 3. — 2500 A luminosity (left) and 2 keV luminosity (right) vs. redshift for the S05 Seyfert 1 (orange upward-pointing triangles), SDSS (Wue circles), and 
high-redshift (green boxes) samples shown with the COMBO-17 (red diamonds), BQS (vioiet downward-pointing triangles), and optically selected, high-redshift 
AGNs (green stars). X-ray upper limits are represented using open symbols (left) or arrows (right). The '2500 A " z SUDSam P le 01 187 sources used in §3.6 is 
denoted by the dotted-line box (Jeff). The properties of the known X-ray-selected AGNs at 2 > 4 (black crosses) are also shown; these sources are not included 
in our formal analyses due to their different selection criteria. 



3.1. Statistical Tools 

While our extended coverage of the I — z plane helps to 
break the redshift bias introduced by using luminosities in- 
stead of fluxes in our analysis, a correlation between I and 
z can still be clearly seen in Figure 3. To take into account 
this existing correlation we use a partial-correlation analy- 
sis method that is designed to measure the correlation be- 
tween two variables, controlling for the effects of one or 
more additional variables. A partial-correlation method that 
properly handles censored data was developed by Akritas & 
Siebert (1996). This method builds upon Kendall's rank- 
correlation coefficient (Kendall 1938) and Kendall's partial 
rank-correlation coefficient (Kendall 1970) to include cen- 
sored data. To measure the partial correlations we use the 
FORTRAN program CENS.TAU, available from the Penn 
State Center for Astrostatistics 9 , which uses the methodol- 
ogy presented in Akritas & Siebert (1996). We present the 
results of our partial-correlation analysis in Table 4. In this 
table we also present the significance of the correlations with 
the potentially obscured (T < 1.6, see § 2.2.1) COMBO-17 
AGNs removed. With the exception of a slight decrease in the 
significance in the Z2kev- '2500 A correlation, there is no sub- 
stantial change in the results when these sources are removed. 
Therefore, we use the full sample in our subsequent analyses. 

To derive the linear-regression parameters for correlations, 
we use the Astronomy SURVival Analysis software package 
(ASURV Rev 1.2; Isobe et al. 1990; Lavalley et al. 1992), 
which implements the bivariate data-analysis methods pre- 
sented in Isobe, Feigelson, & Nelson (1986). ASURV also 
properly handles censored data using the survival-analysis 
methods presented in Feigelson & Nelson (1985) and Isobe 
et al. (1986). We performed linear regressions using both 

9 Available at http://www.astrostatistics.psu.edu/ 
statcodes/cens_tau 



the fully parametric EM (estimate and maximize) regression 
algorithm (Dempster, Laird, & Rubin 1977) and the semi- 
parametric Buckley-James regression algorithm (Buckley & 
James 1979). The EM Regression algorithm uses an iterative 
least squares (ILS) method that reduces to the traditional ordi- 
nary least squares (OLS) when no censored data are present. 
We present the parameters provided by the EM regression al- 
gorithm below, but in all cases the results from the Buckley- 
James regression algorithm agreed within the errors. 

3.2. l 2 kev versus l 250 oA 

Early X-ray studies of AGNs revealed a correlation between 
X-ray and UV monochromatic luminosities, typically found 
to be lx oc Z(j V , where /? ~ 0.7 — 0.8 (see the references in 
§1). The observational evidence that (3 is not unity was un- 
expected and led to suggestions that this finding was not an 
intrinsic property of AGNs, but rather an observational bias 
introduced as a result of the larger scatter of l\jv compared 
to l x (e.g., Chanan 1983; Yuan et al. 1998). La Franca et al. 
(1995) found an lx — luv correlation consistent with (3=1 us- 
ing a generalized orthogonal regression procedure that takes 
into account measurement errors on both variables and intrin- 
sic scatter. However, this method does not properly take into 
account the effects of censored data, which typically affect 
optically-selected AGN samples. Currently there is no sta- 
tistical analysis method that properly handles data with both 
censoring and measurement errors (M. Akritas 2005, private 
communication). 

For our optically-selected AGN sample we find a 
highly significant {15.3a) correlation between the observed 
monochromatic X-ray and UV luminosities. The EM linear- 
regression algorithm in ASURV that we use to calculate the 
linear-regression parameters is based on the traditional ordi- 
nary least-squares method, OLS(F | X), which minimizes 
the residuals of the dependent variable, Y (the EM method 
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10 aB 10 29 1 30 1 31 10 
Z(2500 A) [erg s" 1 Hz" 1 ] 



FIG. 4. — (top) Rest-frame 2 keV monochromatic luminosity versus rest-frame 2500 A monochromatic luminosity. The symmetric, best-fit ^2keV~ ^2500 A 
relationship, given by Equation (lc), is denoted by a solid, black line. Equations (la) and (lb) are denoted by solid, gray lines in order of increasing l 2 keV at 
'2500 A = 10 32 er S s_1 Hz -1 . For comparison, the best-fit line derived by S05 (dashed line) and a/3 = 1 relation (dotted line), normalized to best fit the data, 
are shown. The residuals for Equation (la) and the /3 = 1 relation are shown in the two lower panels, respectively. The residuals for Equation (lb) and the 
fi = 1 relation are shown in the two panels on the right, respectively. The overlaid error bars denote the mean and the 3<r standard deviation of the mean of the 
residuals calculated for each Alog(Z 25 QQ x ) = 1 bin. Symbols are defined as in Figure 3, although all symbols are plotted as open to minimize symbol crowding. 
Limits are denoted with arrows. We find a highly significant (15.3a") correlation between /2keV !m ^ 1 ^2500 A ( conn ' omn g f° r redshift) with a best-fit slope of 
f3 = 0.73 ± 0.1, calculated from the ILS bisector. 



assumes the residuals are Gaussian distributed). Since both 
I2 keV and l 2500 A are observed, neither can be truly called the 
"independent" or "dependent" variable. Indeed, a different 
result can be obtained if the residuals of the independent vari- 
able are instead minimized [i.e., OLS(X | Y)]. Rather than 
use the traditional OLS(F | X), we choose a method that 
instead provides a symmetric fit to the data. Following the ar- 
guments in § 5 of Isobe et al. (1990), we use the OLS bisector, 
which simply bisects the two lines from the OLS(F | X) and 
OLS(X I y),orinourcaseILS(F | X)andILS(X | Y) from 
the EM Algorithm in ASURV. 

We perform linear regressions with ASURV on the full 
sample of 333 AGNs and find the relation between ?2keV and 



2500 A 



to be 



log(Z 2 kev) = (0.64±0.02)log(Z 2500A ) + (6.87±0.63) (la) 



using ILS(F | X) (i.e., treating hh-w as the dependent vari- 
able) and 

log(Z 2ke v) - (0.82 ± 0.02) log(7 250oA ) + (1.74 ±0.75) (lb) 

using ILS(X | Y) (i.e., treating l 2500 a as the dependent vari- 
able). Note that the slope of the I2 keV- '2500 A relation does de- 
pend on which luminosity is used as the dependent variable, 
but both cases are inconsistent with (3 = 1. We use the equa- 
tions given in Table 1 of Isobe et al. (1990) to calculate the 
bisector of the two regression lines. We find the ILS bisector 
to be 

log(/ 2 kev) = (0.72±0.01)log(Z 250oA ) + (4.53±0.69) (lc) 

In Figure 4, we show I2 k e v versus Z 2500 A for the full AGN 
sample (top panel). We show the symmetric, best-fit Z2kev- 



2500 , 



relationship, given by Equation (lc), as a solid, black 
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FIG. 5. — (top) «ox vs. rest-frame 2500 A monochromatic luminosity. The 
best-fit line from Equation (2) is shown (solid line). For comparison, the best 
fit line derived by S05 (dotted line) is also shown, (bottom) Residuals from 
the fit shown in the top panel. The overlaid error bars denote the mean and 
the 3a standard deviation of the mean of the residuals calculated for each 
Alog(/2 5 QQ^) = 1 bin. Symbols are defined as in Figure 3, although all 
symbols are plotted as open to minimize symbol crowding. Limits are de- 
noted with arrows. We find a highly significant (13.6cr) correlation between 
ce x and '2500 A ( controlnn g f° r redshift). 

line. Equations (la) and (lb) are denoted by solid, gray lines 
in order of increasing /2kev at ^500 A = l^ 32 er § s_1 Hz -1 . 
For comparison purposes we show the /2kev- '2500 A relation 
found by S05 (dashed line), along with a = 1 relation with a 
normalization chosen to minimize the /2kev residuals (dotted 
line). The residuals for Equation (la) and the fit assuming 
= 1 are given in the bottom panels, with the mean and the 
(3<t) standard deviation of the mean calculated for sources in 
each Alog(Z 2 5oo a) = 1 bin denoted with large error bars. The 
residuals for Equation (lb) and the fit assuming = 1 are 
given in the right-hand panels using the same symbols. From 
the systematic 0=1 residuals present along both axes, it is 
apparent that a = 1 relation provides an unsatisfactory fit 
to the data, independent of the luminosity considered as the 
dependent variable. 

3.3. ol ox versus l 250 oA 
Given the highly significant correlation between ^keV and 
^2500 A found f° r our AGN sample, we investigate how ot ox 
changes with respect to / 2 50oA' ^2kev, and redshift. For our 
optically-selected AGN sample we confirm a highly signif- 
icant (13.6er) anti-correlation between a ox and Z 2 50oA when 
controlling for the effects of redshift. The best-fit parameters 
for the a ox - l 250Q A relation are 

a ox = (-0.137 ±0.008) log(Z 2500 A ) + (2.638 ± 0.240) (2) 

In Figure 5, we show a ox versus l 2500 \ f° r our AGN 
sample. We show the best-fit linear regression found for this 
sample, given in Equation (2), as a solid line. The residuals 
for the fit are given in the bottom panel, with the mean and 
standard deviation of the mean calculated for sources in each 
Alog(Z 2500 a) = 1 bin denoted with 3er error bars. For com- 
parison purposes we show the a ox - l 250Q ^ relation found by 
S05 (dotted line). 



FIG. 6. — (top) a ox vs. rest-frame 2 keV monochromatic luminosity. The 
best-fit line from Equation (3) is shown (solid line), (bottom) Residuals from 
the fit shown in the top panel. The overlaid error bars denote the mean and 
the 3a standard deviation of the mean of the residuals calculated for each 
Alog(Z2kev) = 1 bin. Symbols are defined as in Figure 3, although all 
symbols are plotted as open to minimize symbol crowding. Limits are de- 
noted with arrows. We find a lower, but still significant (3. Oct) correlation 
between a ox and ! 2 (controlling for redshift). 

3.4. oc ox versus l 2keV 

The correlation between a ox and l 2 k e v has been previously 
examined (Green et al. 1995), but the AGN samples used had 
low X-ray detection fractions and included both radio-loud 
AGNs and BALQSOs, both of which can obfuscate the in- 
trinsic X-ray emission of interest here (see § 1). In addition, 
Green et al. (1995) did not check if the correlation between 
a ox and Z2kev was the result of a true correlation, or a by- 
product of the known correlation between l 2 k e v and redshift. 

The high X-ray detection fraction of our AGN sample al- 
lows us to examine the correlation between a ox and Z2keV- We 
find a weaker, but still significant (3. Oct) anti-correlation be- 
tween a ox and l 2 k e v when controlling for the effects of red- 
shift and taking into account the double-censoring present in 
this relation (i.e., censoring of both the dependent and inde- 
pendent variables). The linear-regression methods used pre- 
viously to derive the parameters of the correlations (EM and 
Buckley-James regression algorithms) are only strictly valid 
if double-censoring is not present. However, given the high 
X-ray detection fraction of our sample, we continue to use the 
aforementioned EM and Buckley-James regression methods, 
assuming the independent variable (l 2 k e y) is detected. Note 
that our choice of linear-regression techniques does not affect 
the significance of the anti-correlation presented above. 

The best-fit parameters for the a ox - l 2 k e v relation from the 
EM regression method are 

a ox = (-0.077 ± 0.015) log(Z 2k ev) + (0.492 ± 0.387) (3) 

In Figure 6, we show a ox versus ^2kev for our full AGN 
sample; the best-fit linear regression found for this sample, 
given in Equation (3), is shown by a solid line. The residuals 
for the fit are given in the bottom panel, with the mean and 
standard deviation of the mean calculated for sources in each 
Alog(Z2kev) = 1 bin denoted with 3cr error bars. 

The clear trend of higher (i.e., less-negative) a ox values 
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Kendall's T m = -0.031 (1.3a) 




ot„(z) = - (0.079 ±0.006)«z-(1.395±0.014) 




1 2 3 4 5 6 

REDSHIFT 

FIG. 7. — (top) a x vs. redshift. The best-fit line from Equation (4) is 
shown (solid line). For comparison, the best-fit line derived by S05 (dotted 
line) is also shown, (bottom) Residuals from the fit shown in the top panel. 
The overlaid error bars denote the mean and the 3cr standard deviation of the 
mean of the residuals calculated for each Az = 1 bin. Symbols are defined 
as in Figure 3, although all symbols are plotted as open to minimize symbol 
crowding. Limits are denoted with arrows. We find no significant (1.3er) 
correlation between a ox and redshift (controlling for '2500 A-*' 

with decreasing Z 2500 A seen in Fig 5 is not as apparent in the 
&ox- h kev plot. This could be due, in part, to the smaller range 
in X-ray luminosity covered by the AGNs (Alog(Z2kev) = 4) 
compared to the UV luminosity range (Alog(Z 2500 A ) = 5). In 
addition, the scatter in the residuals for the a ox - \ i^v relation 
is larger than that seen for i 



'2500 i 



3.5. a c 



versus redshift 

To test for a possible redshift dependence of a ox we mea- 
sure the partial correlation of a ox and z, controlling for l 25O0 x- 
We do not find a significant correlation between a ox and red- 
shift (1.3a), consistent with previous a ox studies. The best-fit 
parameters for the a ox — z relation are 

a ox = (-0.079 ± 0.006) z - (1.395 ± 0.014) (4) 

In Figure 7, we show a ox versus z for our full AGN sample. 
We show the best-fit linear regression found for this sample, 
given in Equation (4), as a solid line. The apparent correlation 
between a ox and redshift is simply an artifact of the Z 2500 A - z 
correlation. The residuals for the fit are given in the bottom 
panel, with the mean and standard deviation of the mean cal- 
culated for sources in each Az = 1 bin denoted with 3er error 
bars. For comparison purposes we show the a ox — z relation 
found by S05 (dotted line). 

In Figure 8, we show the a ox — aoxfcooA) residuals as a 
function of redshift (top panel) and the a ox — a ox (z) residuals 
as a function of Z 2 sooA (bottom panel). It is clear from Fig- 
ure 8 that there is a strong luminosity dependence in the a ox — 
oiox(z) residuals, while the a ox — OoxfeooA) residuals show 
no strong redshift dependence. Using ASURV to find a ox as 
a function of both l 2500 x and redshift, we find the following 
relation. 



: (-0.126 ±0.013) log(/ 250oA ) 
-(0.010 ± 0.009) z + (2.311 ± 0.372) 



Equation (5) shows that the coefficient of z is statistically con- 
sistent with zero (l.ler), which agrees with the findings of S05 
and is consistent with a ox being independent of redshift. 

We use the (\ ox — ^ox(^250oa) residuals to constrain the 
maximum possible residual dependence of a ox on redshift, 
which constrains the maximum amount the ratio of the UV-to- 
X-ray luminosity can evolve with redshift. We find the best- 
fit parameters of the residuals to be <a x — ^oxfeooA) > = 
(-0.004 ± 0.006) z - (0.004 ± 0.012). The maximum abso- 
lute lcr value for the slope of the residuals is 0.010 z, which 
at z = 5 is equal to 0.05. The ratio between the UV and 
X-ray flux is, by definition, r = /„(2500 k)/f v (2 keV) = 
^q2.606q O x Differentiating with respect to a ox , we get Sr/r — 
2.606 ln(10) Sa ox ~ 6Sa ox ~ 0.30, for Sa ox = 0.05. Thus 
the the ratio of UV-to-X-ray flux has not changed by more 
than 30%(1<t) over the redshift range z = — 5. 

To examine if a ox is independent of redshift over the range 
of luminosities observed, we divided the AGNs equally into 
three luminosity bins, with each bin containing 111 sources. 
We plot the a ox — aoxfcoo a) residuals as a function of red- 
shift for each UV luminosity bin in Figure 9. It is apparent 
that the COMBO- 17 sample and the Seyfert 1 sample of S05 
comprise most of the lower-luminosity AGNs (log(^ 2 5oo a) < 
29.812) . The number of COMBO-17 AGNs decreases in the 
moderate UV luminosity bin where the BQS and SDSS AGNs 
become prevalent (27.812 < log(Z 250oA ) < 30.795). At the 
highest luminosities, the z > 4 AGNs and the SDSS sam- 
ple dominate (log(Z 2500j 4) > 30.795). None of the residuals 
shows any indications of having a systematic offset, showing 
that AGNs spanning five decades in UV luminosity show little 
a ox evolution with redshift. 

3.6. Comparing ol ox - z and I2 kev - Z relations 

Recent discussions have questioned the validity of measur- 
ing the redshift evolution of a ox controlling for the effects of 
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FIG. 8. — aox residuals as a function of redshift (top panel) and ' 25 qq A 
(bottom panel). The overlaid error bars denote the mean and the 3<r standard 
deviation of the mean of the residuals calculated for each Az = 1 bin (top 
panel) or Alog(Z 25 QQ^) = 1 bin (bottom panel). Symbols are defined as 
in Figure 3, although all symbols are plotted as open to minimize symbol 
crowding. Limits are denoted with arrows. The systematic residuals in the 
lower plot indicate that a ox cannot be dependent on redshift alone. 
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/2500 A' since ^500 A an d a ox are, by construction, correlated 
(B. Kelly 2005, personal communication). It has been sug- 
gested that measuring the redshift evolution of ^kev, control- 
ling for Z 2 5oo A' ls a better measure of the evolution of a ox since 
this relation does not include the aforementioned a Q x - ^2500 A 
correlation. However, it is not clear that the redshift evolution 
of h keV and a ox are necessarily related, even when controlling 
for the effects of Z 2500 1 (e.g., M. Akritas 2005, personal com- 
munication). We examine the redshift evolution of Z 2 k e v using 
our sample and find a highly significant partial correlation be- 
tween ^2kev - z, controlling for ^500 A ( T i2, 3 = 0.21, it = 8.7). 
If the I2 kev - z partial correlation is indeed a proxy for a ox - z, 
then this result conflicts with our finding in §3.5 that a ox does 
not evolve with redshift. We examine this apparent paradox 
further using both Monte Carlo simulations and a subset of 
our full data sample that better covers the l 2500 a - z plane. 

To determine if the I2 kev - z correlation is real or fake, we 
ran a set of simulations designed to assess the strength of par- 
tial correlations expected for samples with similar Z 2 5oo A " 
z distributions. We started by randomly assigning an / 2 5ooA 
value in the observed range (29.0 <^50oA< 32.5) to each 
AGN in our sample. We then computed l 2 k e v f° r each source 
assuming the ILS bisector relation in Equation lc, including 
a random scatter comparable to that observed. We compute 
the I2 kev - z partial correlation (controlling for Z 2500 ^) for ten 
random samples to obtain an average value for the strength 
and significance of the partial correlation. In these realiza- 
tions there is no l 2500 a - z dependence, by construction. We 
find average correlation strength of 7-12,3 = 0.06 ± 0.02 with 
a = 2.2 ± 0.7, significantly lower than for our actual sam- 
ple. However, if we instead constrain the random / 2 5 00 a t0 



follow the observed ? 25 o a - z dependence, the partial tau 
method finds a more significant, but false, partial correlation 
of Ti2 :3 = 0.17 ± 0.01 with a = 7.0 ± 0.5. These sim- 
ulations suggest that the strength of the correlation between 



h kev and redshift is an artifact of the strong I 
in flux-limited samples. While the partial-correlation meth- 
ods of Kendall's tau are designed to remove the contribution 
of the third variable (in this case / 2 50oa) fr° m the measure- 
ment of the correlation between the first and second variables 
(h kev and z), when the correlations between the controlling 
variable and the first and/or second variable are very strong, 
apparently significant, but false, partial correlations can arise. 

To test further the I2 k e v - z relation, we limit our analysis 
to sources with 29.3 <'2500A< 31 and z < 3, consisting of 
187 sources (see the subsample box in Figure 3). This sub- 
sample more completely fills the l 2500 A ' z plane and thus is 
less affected by the strong ^500 A " 2 correlation present in 
flux-limited samples. Using this subsample, we find ^2kev - z 
partial-correlation results of 7-12,3 = 0.11 with a = 2.8, again 
showing that the strength of the relation diminishes greatly 
when the flux-limited nature of the sample is reduced or elim- 
inated. For comparison, the a QX - ^500 A partial correlation 
decreases, as one would expect from the smaller sample size 
and reduced /2500A and z range, but it remains significant 
( r i2,3 = —0.29, a = 7.0). From these tests it is apparent 
that the redshift evolution of I2 k e v is strongly influenced by 
the strong correlation between l 2500 a and z and thus should 
not be used as a substitute for the a QX - z relation. 

3.7. Investigating non-linear a ox correlations 

From the previous three sections it is clear that a ox is most 
significantly correlated with ^500 A' but does the slope of this 
correlation remain the same over the entire range of UV lu- 
minosities, or is the a QX - / 2 5oo A relation non-linear? From 
Figure 4 it appears possible that the Z2kev- ^2500 A slope may 
become steeper at lower luminosities. It is difficult to tell from 
the plot, however, as lower-luminosity sources typically have 
larger measurement errors. Non-linearity may also be hinted 
at in the residual plot in Figure 5. While all of the 3c error 
bars for the mean residuals in each luminosity bin are consis- 
tent with zero, there is a slight curve seen in the mean residu- 
als, peaking around 10 30 erg s _1 Hz -1 . 

We investigate this potential change in the slope of the a ox - 
'2500 A relation by dividing the sample in half using the me- 
dian UV luminosity (log(Z 2500 ^) = 30.352) and calculating 
the slope of the best-fit line for both the high-luminosity and 
low-luminosity halves of the sample. We find the coefficient 
for the low-luminosity half to be flatter (-0.098 ±0.019) than 
for the high-luminosity half (-0.152 ± 0.019). We use Stu- 
dent's t-test to determine if these two slopes are drawn from 
the same distribution. We find only a 5% probability that the 
best-fit slopes for the high and low-luminosity halves of the 
distribution are drawn from the same parent distribution, sug- 
gesting that the slope of the a ox - Z 2 5oo A relation may be Z 2500 ^ 
dependent. 

We use this same method to determine if the ol Q x~ ^2500 A 
slope is dependent on redshift. We divide the full sample into 
low-redshift (z < 2) and high-redshift (z > 2) subsamples, 
containing 230 and 103 sources, respectively. We do not find a 
significant difference in the best-fit slopes for the low-redshift 
(-0.125 ± 0.012) and high-redshift (-0.147 ± 0.020) sub- 
samples. Using Student's t-test, we find a 37% chance that 
the slopes of the low-redshift and high-redshift subsamples 
are drawn from the same parent distribution. 
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4. DISCUSSION AND CONCLUSIONS 
4.1. Summary of AGN Sample 

In this paper, we examine the correlations between ^500 A' 
hkev, &ox, and redshift for optically-selected AGNs. We 
extend the coverage of the luminosity-redshift plane rela- 
tive to the S05 sample by adding 52 lower-luminosity AGNs 
discovered by the COMBO- 17 survey, 46 high-luminosity, 
low-redshift (z < 0.5) AGNs from the BQS, and 19 high- 
redshift (z > 4) AGNs from recently published, targeted 
X-ray/optical AGN studies. Radio-loud AGNs and BALQ- 
SOs are excluded from our sample whenever possible, and 
the effects of host-galaxy contamination are largely removed 
for the COMBO- 17 AGNs using high-resolution ACS obser- 
vations. These additional AGNs bring the total sample to 333 
optically-selected AGNs, 293 (88%) of which have detected 
X-ray counterparts. 

In Table 5 we present the statistical properties of the full 
sample, broken into Alog(/ 2500 ^) = 1 bins. The Kaplan- 
Meier estimator within ASURV was used to calculate the sta- 
tistical properties of the sources in each bin. The redshift 
range, mean, RMS error, median, and both 25th and 75th per- 
centile values are given, along with the value of Z 2 5ooA an d 
a ox calculated from Equations (lc) and (2), respectively, us- 
ing the mean / 2500 A va l ue given in the appropriate bin. This 
table also gives both the expected value of I2 k e v and a ox for 
each Alog(? 250oA ) = 1 bin as well as the spread around the 
expectation values. From the table, it appears that the spread 
of h kev decreases with increasing Z 2500 A f° r tne ^ ms w hh sta " 
tistically significant numbers of sources. There is no clear 
'2500 A dependence on the spread of a ox . This table is useful 
in identifying AGNs that may be X-ray weak or assessing the 
excess associated with jet-linked X-ray emission. 

4.2. Results of Partial Correlation Analyses 

To measure correlations in our data we employ partial- 
correlation methods that deal with censored data. We confirm 
that the monochromatic 2500 A and 2 keV luminosities are 
highly correlated (15.3ct), considering the effects of redshift, 
and follow the relation I2 kev °c /^ ooA , where = 0.73 ±0.01 
(from the ILS bisector), and not /? = 1 as found in some pre- 
vious studies. 

We investigate how the ratio of the X-ray and UV lumi- 
nosities, a ox , is correlated with Z 250 oA' ^2kev, and redshift for 
our AGN sample. When taking into account the effects of 
redshift, we confirm that the a ox - Z 250oA anti-correlation is 
highly significant (13.6c) but find a much smaller, but still 
significant, anti-correlation for the a ox - hkev relation (3. Oct). 
The discrepancy in the significance of these two correlations 
is due, in part, to the slope of the h kev- '2500 A relation be- 
ing less than unity. As can be seen in Figure 4, the range of 
possible values of I2 kev for a given value of Z 2500 A is smaller 
(and hence the range of a ox values is smaller) than the range 
of possible values of £ 250oA for a given value of ^kev- In ad- 
dition, some AGNs are more variable in the X-ray than they 
are in the optical/UV, with long-timescale X-ray observations 
of Seyfert Is showing variability > 100% for some sources 
(Uttley & McHardy 2004, and references therein), which in- 
creases the possible ranges of I2 kev even more dramatically 
for a given value of l 25QQ \- 

Our analysis did not reveal any significant (1.3ct) evolu- 
tion of a ox with redshift (when considering the effects of 
'2500 a)' which is consistent with the majority of previous 
a ox studies. However, including the COMBO- 17 AGNs in 



TABLE 4 
Partial Correlation Results. 



Relation 


Controlling 
Variable 


Kendall's 

T 12,3 


Significance 




Total Sample 


h keV vs - ^2500 
aox vs. Z 250oA 

aox vs. l 2 keV 
a ox vs. 2 


A z 

z 

z 

hsoo A 


0.519 
-0.377 
-0.080 
-0.031 


15.3cr 
13.6<T 
3.0cr 
1.3a 


Removed COMBO-17 AGNs with T < 1.6 


h keV vs - ^2500 
aox vs. / 250oA 
aox vs. l 2 keV 
a ox vs. 2 


A z 

z 

z 

l 25W A 


0.505 
-0.396 
-0.079 
-0.004 


13.7cr 
13.7cr 
3.0cr 
0.2a 



Note. 

T 12 " T 



Kendall's partial-rank correlation is defined as t 12i 3 = 
, where r xy is the Kendall rank-correlation coefficient between 



[(WsK 1 -^)] 1 ^ 

data vectors x and y (Akritas & Siebert 1996). 

our study allows us to constrain, for the first time, the spec- 
tral evolution of moderate-luminosity AGNs out to high red- 
shifts. Since moderate-luminosity AGNs are the numerically 
dominant type of AGN in the Universe, and X-ray surveys 
have found luminosity-dependent density evolution for such 
AGNs, it is important to understand the spectral evolution of 
these types of sources. The lack of cosmic evolution of a ox 
also suggests that the energy generation mechanisms that cre- 
ate AGN emission locally are already in place at z ~ 5 — 6. 
This agrees with AGN studies which find no significant evo- 
lution in the continuum shape of AGNs at high-redshifts from 
radio (e.g., Petric et al. 2003), optical/UV (e.g., Pentericci 
et al. 2003), and X-ray (e.g., Page et al. 2005; Shemmer et al. 
2005) studies. 

4.3. Comparison with X-ray-Selected, z > 4 AGNs 

In §3, we limited our analyses and discussion to the 333 
optically-selected AGNs in our sample. Here we determine 
if the relations we derived for Z 2500 A , ^kev, aox, and redshift 
also hold for the eight published X-ray-selected, radio-quiet, 
z > 4 AGNs (see §2.2.3). While subject to different selec- 
tion biases than our main AGN sample, these X-ray-selected 
AGNs substantially extend the coverage of the luminosity - 
redshift plane (see crosses in Figure 3a), and thus they can be 
useful in determining if the aforementioned relations are also 
valid for high-redshift, low-luminosity AGNs. 

From Figures 4, 5, and 6 it is apparent that the X-ray- 
selected AGNs, while considerably fainter in l 2500 A than the 
optically-selected AGNs in our sample at similar redshifts, 
have a ox values that do not significantly stand out from the 
a ox values of the optically-selected AGNs. While the X-ray- 
selected AGNs do populate the high (less-negative) portion 
of the a x distributions in Figures 5 and 6, this is not unex- 
pected because X-ray selection is biased toward sources that 
are X-ray bright. In other words, for a group of sources at a 
given ^ 25 ooA with the observed spread in a ox , X-ray surveys 
are biased toward selecting those sources with higher a ox val- 
ues (i.e., the AGNs that are relatively brighter in X-rays). 

In Figure 7, the X-ray-selected sources all lie well above 
the best-fit a ox - z relation derived from our main sample. 
This apparent discrepancy is expected if a ox is indeed depen- 
dent on Z 2 50oA an d not z - The eight X-ray-selected AGNs 
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dox = A [log(/ uv ) - 30.5] + A T [t(z) - 0.5] + A 



This Study 
Strateva etal. (2005) 
AT86 

Wilkes etal. (1994) . 




-0.2 - 



This modification changes Equation (6) of AT86 to 



s 



1 2 

A x 

FIG. 10. — Best-fit values and AS = 4.61 (90%) confidence contours 
for the coefficients A and A T for our sample (black cross, solid black con- 
tour), the sample of S05 (green cross, dotted green contour), and the samples 
examined by Avni & Tananbaum (1986) (open orange circle, dot-dashed or- 
ange contour) and Wilkes et al. (1994) (tilled blue circle, dashed blue con- 
tour). Inset: A magnified view of the contours for our sample. Both the 68% 
(AS = 2.30) and 90% (AS = 4.61) confidence contours (solid black con- 
tours) are shown. For comparison, the 68% and 90% confidence contours of 
the S05 sample are also shown (dotted green contours). 



have considerably smaller / 2500 A values than the optically se- 
lected, z > 4 AGNs, and thus we would predict that they 
would have higher (less negative) values of a ox from Equa- 
tion (2). Indeed, the a ox residuals shown in Figures 8 and 
9 demonstrate that the X-ray-selected AGNs (crosses) agree 
fairly well with our conclusion that a ox is dependent on Z 2500 A 
and shows no significant redshift dependence. To test further 
this conclusion we again ran the partial correlations described 
in §3, this time including the eight X-ray-selected sources. 
We find that the significance of most of the correlations pre- 
sented in Table 4 increase ~ 0.3 — 0.8a. The a ox - z partial 
correlation increases by ~ 0.7<t, but it is still below the thresh- 
old of a significant correlation. In addition, while the X-ray- 
selected AGNs do increase the coverage of the luminosity - 
redshift plane, it is not clear what effect the different selection 
biases of the X-ray-selected AGNs will have on these corre- 
lations. 

4.4. Comparison with Earlier Studies 

The AGN sample presented here enables us to provide the 
tightest constraints to date on the contributions of Z 2500 ^ and 
redshift to a OK . To illustrate our improvement over some past 
studies, we compare our results with those of Avni & Tanan- 
baum (1986, hereafter AT86), Wilkes et al. (1994), and S05. 
Note that the first two authors use cosmological look-back 
time, t(z), in units of the present age of the Universe, and 
not redshift in their calculations. We calculate t(z) for our 
sources and use the method outlined in §§ 3 and 4 of AT86 to 
calculate confidence contours. We rename A z from AT86 to 
A T since, strictly speaking, it is the coefficient of cosmolog- 
ical look-back time which is not a linear function of redshift. 



ox V 6 2500 A' ' 



A [log (I 



2500 



X-ray Loud) = 

^) — 30.5] + At [t (z) — 0.5] + A (6) 



For our sample, we find best-fit values of [A Q) A T , A] = 

[0.140^:^6, -0-009t^4o, L553 ± °- 014 ] for Equation 6. 
In Figure 10, we show the AS* = 4.61 confidence contours 
for our data and the samples presented in AT86, Wilkes et al. 
(1994), and S05. These are equivalent to 90% confidence 
contours taking two parameters to be of interest. The tighter 
constraints provided by our data have a number of sources. 
Our X-ray detection fraction (88%) is higher than either AT86 
(61%) or Wilkes et al. (1994) (64%). Unlike early AGN sam- 
ples used to study a ox , our data represent a relatively homo- 
geneous collection of AGNs where we attempted to remove 
or reduce sources of systematic errors (radio jets, absorption, 
galaxy contribution, etc.). Finally, to our knowledge, our 
AGN sample provides the best coverage of the luminosity - 
redshift plane to date. 

4.5. Comparing X-ray and Optical AGN Luminosity 
Functions 

Understanding the relationship between I2 keV and Z 2500 A m 
AGNs aids in explaining discrepancies between the AGN X- 
ray and optical luminosity functions (XLF and OLF, respec- 
tively). For luminous AGNs, both the OLF (e.g., Boyle et al. 
2000; Croom et al. 2004; Richards et al. 2005, 2006) and 
XLF (e.g., Cowie et al. 2003; Ueda et al. 2003; Barger et al. 
2005; Hasinger et al. 2005; La Franca et al. 2005) are con- 
sistent with pure luminosity evolution (PLE), with the peak 
activity of luminous AGNs occurring at z = 2 — 3. How- 
ever, recent X-ray surveys have found that low-luminosity 
(L 2 -8 kcv 5 10 44 erg s" 1 ) AGNs do not agree with the PLE 
model, but instead appear to undergo luminosity-dependent 
density evolution (LDDE; e.g., Ueda et al. 2003; Steffen et al. 
2003; Barger et al. 2005; Hasinger et al. 2005; La Franca et al. 
2005; Silverman et al. 2005a). These studies found that the 
peak density of low-luminosity AGNs occurs at lower red- 
shifts than for high-luminosity AGNs, an example of "cosmic- 
downsizing". LDDE is not seen in deep optical AGN surveys 
where the OLFs of low-luminosity AGNs are not well con- 
strained and are consistent with both PLE and pure density 
evolution (PDE) models (Wolf et al. 2003). 

Given that AGNs found in deep X-ray-selected samples 
appear to undergo LDDE, does this necessitate LDDE for 
optically-selected samples as well? We can use the LDDE 
equations given in § 5.2 of Ueda et al. (2003, hereafter U03) 
and the relationship between /2keV and ^500 A we ^ n< ^ f° r 
optically-selected AGNs (Equation lc) to examine how the 
LDDE parameters change if the optical luminosity is substi- 
tuted for the X-ray luminosity. Equations (16) and (17) of 
U03 give the evolution factor, e(Lx, z), as a function of both 
the luminosity and redshift, which is unique to LDDE mod- 
els (for PLE and PDE, the evolution factor only depends upon 
redshift). More specifically, the redshift cutoff, z c , included 
in the evolution factor is a function of luminosity in LDDE 
models and is a constant in PLE and PDE models. Apply- 
ing our ^2kev- ^2500 A relation to the redshift cutoff (Equation 
17 of U03), we find that the power-law index, a, becomes 
a' = a ■ (3, where [3 is the slope of the I2 k e v- ^2500 A relation, 

10 In AT86, a ox is defined to be positive. Thus, the coefficients A and A T 
have the opposite sign as the coefficients presented the previous equations. 
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TABLE 5 
Sample Statistics 



log(Z 250() A ) Range 27 - 28 28 - 29 29 - 30 30 - 31 31 - 32 32 - 33 

Number of Sources 6 41 84 121 67 14 

Number of X-ray Detections 5 38 82 99 55 14 

X-ray Detection Fraction 83% 93% 98% 82% 82% 100% 

Redshift Range 0.009 - 0.201 0.009 - 1.264 0.033 - 2.870 0.085 - 4.190 0.472 - 6.280 2.736 - 4.550 



lo ^2500 A) 



Mean log(« 2500 A ) 27.847 28.604 29.622 30.540 31.552 32.259 

RMS log(/ 250Q A ) 0.081 0.250 0.238 0.286 0.292 0.150 

25th Percentile 27.753 28.508 29.453 30.283 31.264 32.120 

Median log(« 2500 A ) 27.825 28.630 29.625 30.609 31.523 32.200 

75th Percentile 27.890 28.781 29.806 30.791 31.742 32.330 



lo g( ; 2kev) 



Mean log(2 2 kev) 


24.571 


25.149 


25.953 


26.472" 


27.105" 


27.600 


RMS log(Z 2k ev) 


0.318 


0.621 


0.431 


0.418 


0.333 


0.348 


25th Percentile 


b 


24.873 


25.792 


26.296 


26.879 


27.340 


Median log(/2 kev) 


24.463 


25.222 


26.018 


26.462 


27.065 


27.730 


75th Percentile 


24.661 


25.622 


26.221 


26.746 


27.337 


27.836 


Value from Equation (lc) 


24.749 


25.235 


25.889 


26.478 


27.128 


27.582 



"ox 



Mean a x 


-1.253 


-1.322 


-1.408 


-1.568" 


-1.708 


-1.788 


RMS a x 


0.093 


0.192 


0.165 


0.198 


0.146 


0.131 


25th Percentile 


-1.343 


-1.451 


-1.493 


-1.657 


-1.810 


-1.865 


Median aox 


-1.309 


-1.291 


-1.402 


-1.547 


-1.691 


-1.786 


75th Percentile 


-1.230 


-1.184 


-1.289 


-1.455 


-1.609 


-1.717 


Value from Equation (2) 


-1.189 


-1.293 


-1.433 


-1.559 


-1.699 


-1.796 



"ASURV changed first datum in bin from a censored point to a detection — mean estimate is biased. 
b No value is reported due to bias by the upper limit and the small number of sources in this bin. 



when the X-ray luminosity is replaced with the optical/UV 
luminosity. The evolution factor becomes independent of op- 
tical/UV luminosity [i.e., e(Ljjv, z) — ► e(z)] only when the 
X-ray and optical/UV emission in AGNs are independent of 
one another (i.e., (3 = 0). Thus, aside from this special case, 
the LDDE observed for X-ray AGN samples implies LDDE 
for optical AGN samples. 

The j3 — 0.725 slope found for the £ 2 50oA~ ^2kev rela- 
tion (Equation lc) evinces that a group of AGNs will span 
a larger range in UV luminosity than in X-ray luminosity. 
In other words, optical/UV surveys must cover a larger lu- 
minosity range than X-ray surveys to probe the same AGN 
populations. This is apparent when comparing the sky den- 
sities of AGNs from deep X-ray and optical surveys (see 
§ 1). The deep X-ray surveys find about an order of mag- 
nitude more AGNs per square degree than the deepest opti- 
cal surveys. Some of this discrepancy is attributable to dif- 
ferent AGN selection effects between the X-ray and optical 
surveys. Host-galaxy contamination is much more prevalent 
in the optical regime, as is intrinsic obscuration. However, 
if we compare the depth probed by the optical and X-ray 
surveys we see that the optical surveys are missing the low- 
luminosity AGNs that do not fit into the PLE paradigm. Stud- 
ies of the XLF find the densities of AGNs with L 2 -8 kcV = 
10 42 — 10 44 erg s _1 peak at lower redshifts than their more 
luminous cousins, a strong indication of LDDE. This range in 
hard X-ray luminosity corresponds to a monochromatic I2 ke v 
range of log(Z 2 kev)= 23.9 - 25.9, assuming T = 2. Us- 
ing Equation (lc), we find the corresponding monochromatic 

/2500 A ran g e t0 t> e l°g(^2500A) = 26.9 — 29.7. From Figure 3a 
it is apparent that the optically-selected AGN samples do not 
extend to faint enough UV luminosities to detect all of these 



AGNs. Even if we assume /? = 1, we find that the optical 
surveys are still too shallow, but the problem becomes less 
severe. Assuming (3 = 1, the range of hkev covered by the 
X-ray selected AGN surveys corresponds to a monochromatic 
Z 2 5oo A ran g e °f l°g(^250oA) = 27.8 — 29.8, an order of mag- 
nitude higher than we get with /? = 0.73. In light of this 
comparison it is not surprising that LDDE is not required to 
describe the evolution of the OLF, because optical surveys are 
not deep enough to probe the AGN luminosities where LDDE 
becomes apparent in XLF studies. 

4.6. Future Studies 

This study presents the highest significance a ox - 1 2500 A cor- 
relations to date, but significant work remains to be done. The 
significance of the measured correlations could be substan- 
tially improved by eliminating the intrinsic scatter in our data 
due to UV and X-ray variability. This could be accomplished 
with concurrent UV and X-ray observations of AGNs with, 
e.g., XMM-Newton. This study could also be extended by 
including more of the most optically luminous quasars known 
at z ~ 1 — 4. The high luminosities of these quasars make 
them visible at all redshifts, and their inclusion would reduce 
the paucity of these sources at low-to-moderate redshifts in 
our current sample. We could also extend our study to in- 
clude Low Ionization Nuclear Emission-line Regions (LIN- 
ERs). These low-luminosity AGNs could better constrain the 
ctox - '2500 A relation, and the additional 1 25 qq a "leverage" at 
low luminosities would help test the possible non-linear a ox 
- 1 2500 A relation suggested by our data. To measure the AGN 
emission in LINERs it is essential to minimize the contamina- 
tion from the host galaxy. To do this, the high-resolution capa- 
bilities of HST and Chandra are needed to isolate the intrinsic 
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AGN emission within LINERs. Finally, there are currently no 
compelling physical models that explain the origin of the ob- 
served a ox - l 2 5oo A relation. Any theoretical model attempting 
to describe the emission mechanisms within AGNs must be 
able to reproduce this observed relation. There is hope that 
attempts to model the environments around accreting SMBHs 
from first principles using relativistic magneto-hydrodynamic 
simulations (e.g., Hawley et al. 1995; De Villiers & Hawley 
2003; Krolik et al. 2005) will be able to explain the observed 
a ox - 1 25 qo A relation, but to our knowledge these simulations 
have not yet found a physical explanation for the observed a ox 

-^2500A relation - 
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APPENDIX 

A. PROPERTIES OF ADDITIONAL z > 4 AGNS 

We list the properties of the publicly available sources from Kelly et al. (2005) in Table Al. To our knowledge, details on 
these sources are not available in the literature. The Chandra data were reduced using the same methods outlined in Vignali et al. 
(2001). In Table A2, we list the properties of the X-ray-selected AGN presented in § 2.2.3. 



TABLE Al 

Properties of z > 4 Quasars Observed by Chandra 



Object 

(1) 


z 

(2) 


N n 
(3) 


^4-Bl450 

(4) 


/2500 
(5) 


log(i 25 oo) 
(6) 


M B 
(7) 


Exp. Time 
(8) 


Count rate 
(9) 


fsB 
(10) 


/2keV 
(11) 


log(yL„) 2 kcV 
(12) 


log(L H B) 

(13) 


C*ox 

(14) 


R 
(15) 


QSO 0910+564 


4.04 


2.73 


20.7 


3.08 


30.98 


-25.9 


22849 


034 +o.i7 


, ,, + 0.67 
1 - J - i -0.4B 


0.99 


44.18 


44.4 


-1.72 ± 0.09 


< 16.6 


PC 1450+3404 


4.19 


1.26 


21.0 


2.33 


30.89 


-25.7 


14842 


094 +0-32 
"■^-0.25 


, 40+1.21 

-0.93 


2.71 


44.64 


44.9 


1 ci+0-08 

- L " JJ --0.07 


< 20.3 


SDSS 1413+0000 


4.08 


3.12 


19.8 


7.05 


31.35 


-26.8 


11843 


, , 7 +0.40 
■ —0.32 


4 -61 + 125 


3.49 


44.74 


44.9 


-1.65 ± 0.07 


< 7.4 


SDSS 0050-0053 


4.33 


2.69 


19.5 


9.73 


31.53 


-27.3 


12737 


1 i s + 0.39 
1-1 -0.29 


4 24+1.44 


3.38 


44.75 


44.9 


-1.71 ± 0.07 


< 5.6 


SDSS 1444-0123 
SDSS 2357+0043 


4.18 
4.36 


3.92 
3.28 


19.5 
19.8 


9.64 
7.38 


31.50 
31.41 


-27.2 
-27.0 


10001 
12657 


0.40+° o ;f 9 
0.85+°; 3 2 i 


'• 60 -O.76 
J '^'-1.03 


1.24 

2.73 


44.30 
44.67 


44.5 
44.9 


-1 87+ 011 
1 - a ' -0.12 

-1.70 ± 0.08 


< 5.4 

< 7.0 



NOTE. — (1) Source name; (2) Redshift; (3) Galactic column density, from Dickey & Lockman (1990), in units of 10 20 cm -2 ; (4) AB magnitude at rest-frame 1450 A; (5) 
flux density at rest-frame 2500 A, in units of 10 -28 erg cm -2 s _1 Hz -1 ; (6) log of the monochromatic luminosity at rest-frame 2500 A, in units of erg s _1 Hz -1 ; (7) absolute 
B— band magnitude; (8) Chandra exposure times (in seconds) corrected for detector dead time and high-background periods; (9) observed count rate computed in the soft band, in 
units of 10 -3 counts s — 1 ; (10) Galactic absorption-corrected flux in the observed 0.5-2 keV band, in units of 10 _1 ° erg cm -2 s — 1 ; (11) rest-frame 2 keV flux density, in units of 
10 -32 erg cm~ 2 s _1 Hz -1 ; (12) log of the luminosity at rest-frame 2 keV, in units of erg s _1 ; (13) log of the 2-10 keV rest-frame luminosity, corrected for the effect of Galactic 
absorption, in units of erg s _1 ; (14) optic al-to-X-ray spectral index; errors have been computed following the "'numerical method" described in § 1.7.3 of Lyons (1991); both the 
statistical uncertainties on the X-ray count rates and the effects of the observed ranges of the X-ray and optical continuum shapes have been taken into account; (15) radio loudness 
parameter, defined as R = ghz// 440Q ^ (rest frame). The 5 GHz flux density is computed from the 1.4 GHz flux density (from FIRST) assuming a radio power-law slope of 
a — —0.8, with f v oc v a ; 



TABLE A2 

Properties of X-ray-Selected, z > 4 AGNs 



Object 


z 




JV H 


^1-Bl450 


/2500 


log(L 2 50o) 


/SB 




/2koV 


log(L2kcv) c* ox Reference 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 




(7) 




(8) 


(9) (10) (11) 


CXOCY J033716.7-050153 


4.61 


4.82 


23.8 


-28.84 


29.74 


1.8 


-31.82 




26.76 


-1.14 


Treister et al. (2004) 


CLASXS 11032414.2+572227 


5.40 


0.73 


24.3 


-29.04 


29.65 


0.85 


-32.09 




26.60 


-1.17 


Steffen et al. (2004); Yang et al. (2004) 


RX 11052+5719 


4.45 


0.56 


22.6 


-28.36 


30.20 


2.3 


-31.73 




26.83 


-1.29 


Schneider et al. (1998) 


CXOMP 1105655. 1-034322 


4.05 


3.55 


22.3 


-28.24 


30.26 


5.2 


-31.37 




27.13 


-1.20 


Silverman et al. (2005b) 


CXOHDFN 1123647.9+620941 


5.19 


1.48 


23.5 


-28.72 


29.94 


0.29 


-32.57 




26.09 


-1.48 


Barger et al. (2002); Vignali et al. (2002) 


CXOHDFN 1123719.0+621025 


4.14 


1.45 


25.2 


-29.40 


29.11 


0.26 


-32.70 




25.81 


-1.26 


Barger et al. (2002); Vignali et al. (2002) 


CXOCY J125304.0-090737 


4.18 


2.96 


23.0 


-28.52 


30.00 


1.7 


-31.88 




26.64 


-1.28 


Castander et al. (2003) 


CXOMP 12 1 3945.0- 234655 


4.93 


3.55 


21.6 


-27.96 


30.67 


1.8 


-31.80 




26.83 


-1.47 


Silverman et al. (2002) 



NOTE. — (1) Source name (in order of increasing Right Ascension); (2) Redshift; (3) Galactic column density, from Dickey & Lockman (1990), in units of 10 20 cm -2 ; (4) AB 
magnitude at rest-frame 1450 A; (5) log of the flux density at rest-frame 2500 A, in units of erg cm - 2 s — 1 Hz - 1 ; (6) log of the monochromatic luminosity at rest-frame 2500 A, in 
units of erg s~ 1 Hz - 1 ; (7) Galactic absorption-corrected flux in the observed 0.5-2 keV band, in units of 10~ 15 erg cm -2 s~ 1 ; (8) log of the rest-frame 2 keV flux density, in units of 
erg cm -2 s _1 Hz -1 ; (9) log of the monochromatic luminosity at rest-frame 2 keV, in units of erg s _1 Hz -1 ; (10) optic al-to-X-ray spectral index; (11) References for X-ray-selected 
AGN. 



